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■ We derive an expression for the drag rate (i.e., interlayer momentum trans- 
fer rate) for carriers in two coupled two-dimensional gases to lowest nonvan- 
ishing order in the screened interlayer electron-electron interaction, valid for 
arbitrary intralayer scattering mechanisms, using the Boltzmann transport 

£h ' equation. We calculate the drag rate for experimentally relevant parameters, 

■ and show that for moderately high temperatures (T > 0.2Tp, where Tp is the 
Fermi temperature) the dynamical screening of the interlayer results in a large 

q \ enhancement of the drag rate due to the presence of coupled plasmon modes. 

This plasmon enhancement causes the scaled drag rate to have a peak (i) as 
a function of temperature at T ~ 0.5Tp, and (ii) as a function of the ratio of 

■ densities of the carriers in the two layers when their Fermi velocities are equal. 
We also show that the drag rate can be significantly affected by the intralayer 
scattering mechanisms; in particular, the drag rate changes approximately by 
a factor of 2 when the dopant layer modulation doped structures are moved 
in from 400 A to 100 A. 
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I. INTRODUCTION 



Coupled quantum wells fabricated by epitaxical growth make up interesting systems for 
studying electron-electron interactions in the low- dimensional systems both experimentally 
and theoretically. In coupled wells placed within tunneling distance, one can for instance 
study Coulomb gaps in high magnetic fields!] For wells separated so that tunneling can 
be disregarded, the effect of the mutual polarization may become important even at zero 
magnetic field, and it has been suggested that the interlayer correlations can drive a Wigner 
crystallizationB when the coupled plasmon mode goes soft. This is however only relevant for 
low density systems and has so far not been observed experimentally. At higher densities, 
coupled collective modes in multilayer systems have been seen in inelastic light scattering 
experimentsSB In this paper, we discuss another probe of interlayer interactions and the pos- 
sibility of observing the coupled plasmon modes in bilayer systems, namely by the so-called 
Coulomb drag effect ,!'! where a current in one layer drives a current in the other layer due to 
the momentum loss caused by interlayer electron-electron scattering events. The Coulomb 
drag effect has recently attracted much experiment a]@~0 and theoretical attention.0~lll 

Normally, the effects of electron-electron collisions only have indirect consequences for 
transport properties of single isolated quantum wells, because they conserve momentum. The 
Coulomb drag effect is unique in that it provides an opportunity to directly measure electron- 
electron interaction through a transport measurement where momentum is transferred from 
one layer to the other. Since interlayer interaction well depend strongly on the many- 
body effects of the system it is therefore essential to include screening in the theoretical 
understanding of the measurement s.EH^l This has been realized in much of the theoretical 
work so far. Recently, we have pointed out that at intermediate temperatures on the scale of 
the Fermi temperature the drag effect should in fact be greatly enhanced by "anti-screening" 
due to coupled plasmon modes.E-3 

In this paper, we offer a more detailed study of the predicted plasmon enhancement 
effect in zero magnetic field. We base our calculation on the Random Phase Approximation 
(RPA) which is generally believed to be valid for not too low densities such as e.g. the 
experimental realization fabricated by Gramila et a/.i However, for low density systems 
such as the electron-hole system studied by Sivan et a/.,0 the RPA approximation is not 
able to account for the magnitude of the observed drag effect. This discrepancy between 
theory and experiment could be an indication of break-down of RPA due correlation effects.0 
The Sivan et al. system is however not relevant for the plasmon enhancement effect studied 
here due to the large mismatch between the Fermi velocities,0 as will be explained below. 

The results obtained in this paper are based on the Boltzmann equation. Quantum 
effects neglected by the Boltzmann equation can be shown to be experimentally negligible 
in the structures we consider. 0,0 

We derive a generalized formula for the drag rate which is 
valid for arbitrary intralayer scattering processes, which reduces to the previously obtained 
resultlffi for the case of a constant relaxation time approximation. 

The role of intralayer carrier-carrier interactions has not been studied so far. As men- 
tioned above the intralayer interactions usually only have indirect consequences for the 
transport properties. However, it turns out that the plasmon part of the drag effect is very 
sensitive to the details of the distribution function. We have solved for the distribution func- 
tion including intralayer electron-electron interaction using parameters for GaAs-AlGaAs 
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quantum wells and found that: (1) for high mobility samples where the distant to the 5- 
doped impurities is fairly large the distrubution is well-described by a shifted Fermi-Dirac 
distribution function and which allows us to assume an energy independent momentum 
relaxation rate when calculating the drag rates. This is the situation previously studied 
by several authors.i'EHl^ However, for the case where the dopants is moved closer to the 
two-dimensional electron gas the situation is changed and we predict and enhancement of 
the drag rate because of a larger plasmon contribution in this case. The reason for this 
enhancement is that distribution function becomes skewed toward higher energies and the 
overlap with the plasmon branch consequently becomes larger. 

The paper is organized as follows. In Sec. [TT] we discuss the general formula for the 
drag rate, and in Sec. |T| we study the RPA approximation of the screening for a coupled 
two-layer system. The drag rate for the case of shifted Fermi-Dirac distributions, which is 
applicable when the charged dopants are far from the quantum wells, is studied in Sec. [IV. 
In Sec. |V|, we study the case when distribtions deviates from shifted Fermi-Dirac functions, 
as when the dopants are close to the quantum wells. Finally, a summary is given in Sec. VI . 
Technical details are given in a number of appendices. 



II. BOLTZMAN EQUATION CALCULATION OF THE DRAG RATE 

In a drag experiment, two doped and independently contacted quantum wells are placed 
very close to each other.! We consider the case typical in experiments where one applies 
electric field E\ causing a current of density J\ to flow in layer one. This sets up an induced 
electric field E 2 in layer two, where the current is set to zero. (Throughout this paper we 
denote layer 1 as the driven layer and layer 2 as the dragged layer.) The transresistivity 
tensor pff is then defined as 



P2lJl,a = E 2 ,/3 (1) 



8a 
a=x,y 

By the Onsager relationship, = P21 ■ 

In the case of isotropic parabolic bands, where one has a well-defined effective mass m, 
one can define a drag rate r 2 ~\ through 

- (2) 

where n\ is the carrier density in layer 1. Physically, the drag rate is the net average rate of 
momentum transfered to each particle in layer 2, per unit drift momentum per particle in 
layer 1; i.e., 



1_ = (dp 2 /dt) 
T21 Pi 

where Pi is the momentum per particle in layer i, and the overbar denotes an ensemble aver- 
age. Note that by this definition, the drag rate is not symmetric; i.e., r^ 1 = m 2 n\m{ x n 2 r 2 i, 
which is important when layers 1 and 2 are not identical. We follow the convention of pre- 
vious authors by denoting r 2 \ as rp. 
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The drag rate has previously been derived using a Boltzmann equation approacbMllj as- 
suming a energy-independent intralayer momentum relaxation time, and also using a memory 
functional method with a constant relaxation time approximation. This result was recently 
generalized, using the Kubo formalism, to include energy-dependent impurity scattering 
rates.0 Here we shall further generalize this result to arbitrary intralayer scattering mecha- 
nism within the linear response Boltzmann equation description. 



A. Boltzmann equation for coupled quantum wells 

We define the function ip(k) which is related to the deviation of the distribution function 
from equilibrium by 

Sf(k) = f(k) - f(k) = f(k)[i - f(kM(k) 

- - k ° T PS*) m - (4) 

where f°(k) is the equilibrium Fermi-Dirac distribution function. With this definition, the 
linearized interlayer electron-electron collision term i 



f dki f dq , 

S[i>uihm) = 2j j—^ J j—^ ^(Q^k 1+ q- £ k^ 



f°(k{)f°(k 2 )[l - ^(fci + q)][l - f 2 {k 2 - q)} 
x [^i(fei) + ^ 2 (k 2 ) - i)x{k x + q)- ^ 2 {k 2 - q)} 

6 { £ k 1 + £ k 2 - £ k 1+ q- £ k 2 -q) ■ ( 5 ) 

Here w(q,u) is the probability of a particle scattering with change of momentum and en- 
ergy of hq and tkv, respectively. This rate is usually taken to be the Born approximation 
result w(q,uj) = 2nh~ 1 \Wi 2 (q, w)\ 2 , where Wi 2 (q,u>) is the dynamically screened interlayer 
Coulomb interaction matrix element .0 

We assume that the interlayer interactions are weak, so that these interactions are only 
kept to lowest order. Then, within linear response to an external driving field E\, the coupled 
Boltzmann equations for the system reads 

e 1 E 1 -v(?£±-) = -H 1 [1> 1 ](k 1 ), (6a) 



de 

e 2 E 2 ■ v (^) = S[^ 2 = 0](fe 2 ) - H 2 [V> 2 ] (k 2 ). (6b) 

where ej is the carrier charge in layer i, Hi is the negative of the linearized miralayer 
collision operator.El We have used the assumption of weak interlayer interaction to neglect 
the interlayer electron-electron collision term in Eq. (|6"a| ) and to set ip 2 (k 2 ) = in the 
interlayer collision term S in Eq. (pbD, as these terms are higher order in the interlayer 
interaction than the other terms in their respective equations. For convenience, hereafter 
we shall refer to S[ipi,ip 2 = 0}(k 2 ) simply as S(k 2 ). 

The formal solutions for ip\{k\) and ip 2 (k 2 ) for E fields in the ^-direction are 
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ih(ki) = -e x E x H? 
^ 2 (fc 2 ) = e 2 E 2 H 2 ~ l 



' (dft\ 



(fci) 

(fc 2 )+i7 2 - 1 [S](fc 2 ). 



(7a) 
(7b) 



We assume the electric field in layer 2 is adjusted so that the current in layer 2 is zero, i.e., 



32,x = -2e 2 k B T 



dko 



(^)V*(* 2 )=o, 



(27r) 2 ^V,96 
which implies, by substitution of Eq. ([7b|) into Eq. 



(8) 



2k B Te\E 2 



f dk 2 /<9/ 2 °x -y_ x T / 



2e 2 fc B T 



(fe 2 ) = n 2 e 2 fi t , 2 E 2 
dk 2 r dfl 

V x ,2 ' 



(2tt 



5c 



(9) 



The equivalence in Eq. (|S]) comes from Eq. ( |A5| ) in Appendix |A], and the /ij )2 is the mobility 
of layer 2 in the absence of interlayer coupling. 

Substituting Eq. @ (with ip 2 {k 2 ) = 0) into Eq. (H), and using the identities 



5 (efei + £ k 2 - £ k 1+ q ~ £ k 2 -q) = h y_ oo ^ 5 ( £ fe! " " M%fc 2 " £ k 2 q + M; (10a) 

- /°(e 2 )] = [/°( £2 ) - / ( £l )K( £l - e 2 ), (10b) 

one obtains 

dq 



n 2 e 2 n tj2 E 2 = -Ae 2 %k B T J 
dki 



(27T) 



/oo 
do; w) n B (huj) n B {—huj) 
-oo 



[^(fcj) - Mki + q))[fKki) - A°(fei + g)]*(e fcl - e fel+g " M 



(2tt) 2 
rffc 2 



^(^)^ - / 2 °(fc 2 - «)} 5(e k2 - e k2 + M] }■ (11) 



To write Eq. (|TT|) in a more tractable form, we define the transport relaxation time of 
layer i, Tj(fej), for arbitrary intralayer scattering by 



4>i(ki) 



-eEi ■ Hr 1 



V; 



de 



(ki) 



dEj ■ v(kj) Ti(kj) 
k B T 



(12) 



in analogy with the case of impurity-dominated scattering (see Appendix |A|). Then, from 
Eq. ([7a|) and the definition in Eq. (12), the term in the large square brackets in Eq. ([TT|) 
(i.e., involving integration over ki) is equal to 

/dki 
—± [^i(fei) - Mki + 9)][A°(fei) - /i(fci + q)]6(e ki - s ki+q - hu>) 



e\E\ r dki 



[v x ,i(ki + q)ri(fei + q) - v x ^ 1 (k 1 )T 1 {k 1 )} 



k B T J (2tt) 2 
x [/?(**) - y?(fci + q)]5(e ki - e ki+q - hw). 



(13) 



Note that this quantity is odd with respect to luEi 

Furthermore, using the fact that H 2 X is Hermitian (see Appendix |A|) and Eq. fll2|) , the 
term in Eq. (|11|) in the large curly braces (i.e., involving integration over k 2 ) can be rewritten 

as 

(••jinEq. (11)= {f 2 {k 2 )-f 2 {k 2 - q )}5{e k2 -e k2 _ q + huj) 

x4- 1 K 2 (^)](fc 2 ) 

dk 



l J^k {f"^-f"^-^}^ £ k 2 -^k 2 - q + ^Mk)T2(k) 



k B T J (2vr 



2k B T J (2tt) 2 



x 



S ( £ k' 2 ~ e k' 2 +q ~ hiU ) v xA k 2 + Q)M k 2 + l) 
+6 ^k' 2 ~ £ k' 2+ q + hu)v x ^k' 2 )T 2 (k' 2 )\ . (14) 

For the last equality in Eq. (Ill), we have symmetrized by dividing into two equal parts and 
performing the variable changes: k 2 = k' 2 + q in the first term and k 2 = —k 2 in the second 
term and using inversion symmetry in fc-space. Inserting Eqs. ( |TB| ) and (13) into Eq. (|TTJ), 
and performing a variable change uo — > —uo in the second term in the square brackets of Eq. 
(|H|) (which gives an overall minus sign because of the oddness of Eq. ( |T3D with respect to 
u) gives the following symmetric form for the transresistivity tensor 

£a E 2 h 2 I dq r°° ^ \W 12 {q,u)\ 2 F?{q,u)}F${q,u) 

riiei/Jt^Ei 2ireie 2 n 1 n 2 k B T J (2tt) 2 Jo smh 2 (hu/[2k B T]) 

Here, we have used An B {Tiuj)n B (—huo) = — sinh~ 2 (ftu;/3/2) and we have defined the function 

F°{q, 00) = ^ / j^y 2 [f°(k) - f°(k + q)] 6(e(k) - e(k + q) - hu>) 

x[v a (k + q)T(k + q) -v a (k)T(k)}. (16) 
As required by the Onsager relation, pf? = P21 i n Eq. fllpp. 



B. Drag rate for isotropic parabolic bands 

We concentrate on the case where the bands are isotropic and parabolic, as in the case 
of the electrons or low-energy holes in GaAs, generally the material of choice so far in 
performing drag experiments. In this case, when one has a well-defined mass, the mobility 
can be written in terms of the transport time 

\h = • (17) 

m 

This transport time is the time that enters the conductivity, a = ne 2 T tT /m. One can write 
F a (q,u) = q a Y(q,uj), where 
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Q T tr,i 



rlk 

{/,(£(*)) - fi(e(k + q))} 



(2tt 



x {kn{e{k)) - (k + g)r i (e(fe + g))} 



<J(e(fc) + q) - a) 



(18) 



Hence, using Eq. @ and Eq. (15), the drag-rate for isotropic parabolic bands can be written 

as 



oo 



h f 00 

= T21 = o , rp 2 / ^ ^ / ^ 

|Wi2(g,^)| 2 yi(g,^)y 2 (g^) fl9l 

8wh 2 {hu}/[2k B T\) ' 1 J 

When Ti(k) is independent of k the function V reduces to the usual imaginary part of 
the polarization function, Im[x], and the drag rate formula in Eq. (|19"D reduces the result of 
Refs. [| [TO] and |T2|. We emphasize that results in Eq. ( |7^ ) and Eq. ( |7^ ) ZioZd /or arbitrary 



intralayer scattering mechanism. However, in order to compute the drag rate for a given 
system one must obtain r{k) by first solving for the linear-response distribution function for 
each layer in the absence of interlayer coupling. 

In this paper, we study the drag effect at intermediate temperatures where phonon 
scattering provides an important contribution to the scattering rates and hence one may 
expect that the solution of the Boltmann equation does give a non-constant r,(fc). On 
the other hand, electron-electron interactions tend to pull the distribution function back 
to a shifted Fermi-Dirac distribution. We have performed a numerical calculation of the 
distribution function in GaAs based quantum wells with dopants placed approximately 700 A 



away, as in previous experiments (e.g., in Refs. g and [10]) and found that the distribution 
functions in these cases can, to a very good approximation, be described by a constant 
scattering time. Details are given elsewhereSlj The characteristic time scale for electron- 
electron scattering is given byii r e _ e w {h/Ep) (Ep/kBT) 2 . For GaAs doped at 1.5 x 10 11 
cm -2 at 10 K, r e _ e w 10~ n s, as compared to the impurity and phonon scattering times 
which are typically of order Ti mp> ph ~ 10~ 9 s. The relation r e _ e <C Ti mPiP h explains why the 
shifted Fermi-Dirac distribution function is a good approximation to the actual solution of 
the Boltzmann equation. 

Therefore, in section we shall restrict ourselves to the case where the V-function in 
Eq. (18) can be replaced by the imaginary part of the polarizability|00 



rn 1 



POO POO 

/ dq q 3 dm 
Jo Jo 



8min 2 A;sT7r 2 Jo Jo 
|H / 12 (g,a;)| 2 Imxi(q,^)Imx 2 (q,^) ^ 

However, we show in section that in the case when the impurities in the modulation 
doped samples are moved in closer to the quantum wells, the momentum relaxation times 
r(fc) become significantly dependent on k, and there are fairly large differences between the 
drag rates evaluated by Eq. ( |Hf ) and Eq. fl20p. 
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III. COUPLED PLASMON MODES IN THE RANDOM PHASE 

APPROXIMATION 



To calculate the drag rate, one needs to evaluate the dynamically screened Coulomb 
interaction W(q,uj), which we do using the RPA. The RPA equations for the coupled two 



layer system readi 





VnXi V 12 X2 

V 2 lXl V22X2 




(21) 



where the matrix Vy defines the unscreened Coulomb interaction given by Vy = FijV(q). 
Here, V(q) = 2ne 2 /(eoq) = littFqTF / '(k 2 F q) and Fy are the form factors, for which we use 
the form for square wells (see e.g. Ref. |12]), 

_ 3x + 8n 2 /x _ 327r 2 (l-exp(-x)) 

x 2 + 4tt 2 a; 2 (a; 2 + 47r 2 ) 2 ' 1 ' 
_ 647r 4 sinh 2 (q;/2) ■ 

Fl2 ~ x 2 (x 2 + 4vr 2 ) 2 6 • (22b) 

where x = qL, L is the width of the quantum wells (equal widths are assumed) and d is the 
center-to-center well separation. The solution for the screened interlayer interaction thus 
becomes 

U l2 (q,u) = ^\, (23) 

where 

e(q,u) = [1 - V n (q)xi(q,u)][l - V 22 (q)x2(q^)} 

-[Vi 2 (q)} 2 Xi(q,uj)x 2 (q,uj) (24) 

The collective modes of the coupled electron gas are given by the zeros of the dielectric 
function. There are two such modes, one where the electron densities in the two layers 
oscillate in phase, which we call the optic mode, and one where the oscillations are out of 
phase, referred to as the acoustic mode." At zero temperature U± 2 (q, uj) has two poles on the 
real o;-axis, giving (^-function peaks. At finite temperatures these poles move off the real u> 
axis, which implies that they will gain finite widths. In calculating the plasmon dispersions 
at finite T, we search for the solutions of Re[e(q,uj(q))} = 0; so long as the damping is small 
and the poles are not too far from the real axis, this criterion accurately gives the position 
of the pole. 

While there are simple analytic expressions for Im[x(g, a;)] and Ke[x{q, &)] at T = 0, none 
exist at finite temperatures. Since it is necessary to evaluate the x a ^ finite temperature 
(using the T = expressions for x do not give the required plasmon enhancement, as 
explained later), we have developed an efficient way of calculating the finite-temperature Xi 
which is elucidated in App. 0. 

In order to gain insight into the structure of the collective mode dispersions, we take in 
the following simplifying limit which allows us to obtain an analytic solution of the dispersion 
relation of the the poles at T = 0. 
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A. The thin layer or small-g limit of e(q,uj) 



In the case when qL <C 1, one can to a good approximation set x = qL = in Eq. 



and the form factors reduce simply to F\\ = 1 and F\% = exp(— qd). Then, the dielectric 
function for two identical layers can be written as 

e(q, co) = (1 - e-^ d ) [ X (q, u) - *+(?)] [ X (q, «) - X-(<7)]M?)} 2 , (25) 

where 

X±(ff) = (26a) 

c ±(?) = YlTtF^' (26b) 



J26 



From this it is clear that the dispersion of the two collective modes are given by the conditions 
x(<ly UJ ±(q)) = X±{o)- The zero temperature real part of x is given bj 

q T Fk F ( Q 



V(q) X '(q, w) = ^ (j£ - V 4 " 1 + V a - - 1 J > ( 27 ) 

where we have defined a± = u/vpq ± q/2kp, and used the notation x = x' + We then 
find for the plasmon solutions 



b±Jl + 4/A± 

W± (g) = gw F , (28a) 

fc t = c t -^ r + f, (28b) 

F™F ™F 

a ^SA i+ ^)- (28c) 

The small-g behaviors of the plasmon solutions are 

« y/qqrFVF, (29a) 

in agreement with Ref. Thus the optical mode has the familiar square-root behavior 
while the acoustic mode is linear in q. Because the latter is lower in energy, it gives the 
dominant contribution to the drag rate at small temperatures, as we show later. 

Fig. m shows the dispersion relation of the plasmons as given by Eq. (|28a ), for two different 



values of layer sepation d. Also shown is the part of the q — to plane where particle-hole 
excitations are allowed (i.e., Im[x(g, u)] ^ 0) at zero temperature. 

In the next section, however, we will see that the plasmon contribution to the drag rate 
is in fact given by g-values in the intermediate range, and thus in order to get a quantitative 
theory of the collective mode enhancement, use of the small q expansion is not sufficient, 
and it is necessary to perform a full numerical calculation. 
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IV. THE DRAG RATE FOR THE SHIFTED FERMI-DIRAC DISTRIBUTION 



As mentioned previously, when the dopants are far enough away from the quantum well 
(> 400 A), the distribution functions are to a very good approximation shifted Fermi-Dirac 
functions, and one can use Eq. (p0|) to evaluate r^ 1 . We do so in this section. 



A. Plasmon-pole approximation 

First, we develop a pole approximation for the plasmon contribition to the drag rate and 
compare it with the full numerical solution. For frequencies near the zeros of the real part 
of the dynamical dielectric function we may approximate e(q, uj) = e' + ie" in Eq. (|25|) as 

\e(q,u)\^2V(q)e^ d \i3 ± (q)(u-u ± (q)) + i X "(q,u; ± )\, (30a) 
dx'(q,u 



#fc(ff) 



duo 



(30b) 

UJ=UJ± 



Inserting the approximate dielectric function into Eq. ( |20|) and assuming that the imag- 
inary part of x is small, we approximate the Lorenzian by 5-functions and obtain 



e(q,u) 



5(u-u ± (q)), (31) 



4Im[ X (q,u;)]|/? ± (q)| 
which leads to the following two plasmon contributions to the drag rate 
1 h 2 f«c,± 3 Im[x(q,w±(q))] 



dqq 6 LAV V (32) 

The parameter q c in this expression defines the value of q where the plasmon ceases to 
exist. Operationally, we take g Cj ± to be the wavevector at which the plasmon dispersions 
disappear (i.e., there are no more solutions Re[e(g, u)} = 0) in the finite-T RPA formalism. 
Interestingly, the exponential dependence on the well separation d has dropped out of the 
integrand in Eq. fl3"2"|). However, r± is still (^-dependent, through the ^-dependence of uj±(q), 
(3±(q) and q c . As d is increased, the slope of oo_(q) increases and q c decreases, as evidenced 
in Fig. [I] (at finite temperatures, the acoustic and optic plasmons merge when they come 
close, and the merging point gives q c ). Both these effects tend to decrease the integrand and 
hence the drag rate. On the contrary, the decrease of w + (q) with increasing d tends to make 
t^, 1 bigger, but this is generally a weaker effect. Hence, as expected, r^ 1 decreases with 
increasing well-separation. Obtaining a precise analytic functional dependence of r^ 1 on d 
is difficult, but we have shown numerically that in the plasmon-dominated region r^ 1 ~ d~ a 
where a ~ 3.0 Also note that Eq. (^) breaks down for large d because in this limit, oj + {q) 
and u)-(q) come so close together that the Lorenzians from both these lines overlap. This 
contradicts the initial assumption that the Lorenzians individually can be approximated by 
5-f unctions. 

From Eq. (^) we see that the plasmon contribution is given an integral over, among 
other factors, the imaginary part of x f° r frequencies and wave- vector corresponding to the 
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plasmon dispersions. Therefore a non-zero Im[x] at u±(q), which is outside the T = 
particle- hole continuum, is necessary to obtain any plasmon enhancement effect. Since the 
T = form of Im[x] is always zero at u>±(q), it will never give a plasmon enhancement. One 
must therefore use the finite-temperature form of \ in the evaluation of the drag rate. 

At small temperatures, Im[%] at w±(q) is small because the carriers generally do not 
have sufficient energy to be excited very far above the Fermi surface. However, at inter- 
mediate temperatures (on the scale of the Fermi-temperature) there are enough thermally 
excited particles so that Im[x] at u)±(q) is large and the plasmons actually dominate the 
drag response. 

In Fig. |2] we show the integrand of the Eq. which gives the drag rate at two different 
temperatures (for layers of zero width). There are pronounced peaks are caused by the 
presence of plasmons due to the dynamically screened interlayer interaction. We also show 
the plasmon-pole approximation to the integrand described above, which indicates that this 
approximation is very good for temperatures less than 0.57V f°r the relevant q's which 
contribute significantly to the integral. 

In Fig. |3| we show a contour plot of the integrand of Eq. fl20"|) for T = 0.37V. Note that 
the integrand has significant weight at intermediate q and u values, i.e. up to q ~ .8kp and 
Kuj ~ 2Ep, showing that small u- and g-expansions are generally inadequate for calculating 
T^ 1 at intermediate and high temperatures. 

B. Numerical evaluation: r^ 1 

We evaluated r^ 1 for parameters corresponding to high mobility GaAs samples similar 
to those of ref. |8|. We did this several different ways. First, we numerically integrated Eq. 
(p0|), using the finite-temperature RPA Im[x] obtained through the technique described in 
Appendix [B|. Fig. f| shows the result of this full numerical integration of Eq. (p0|) (solid 
line). We also solved for the temperature dependent plasmon poles numerically and used 
the plasmon-pole approximation described above to obtain plasmon contribution to the drag 
rate. One can see that the approximate curves reproduce to a certain degree the enhance- 
ment seen in the full solution. At low temperatures the upturn in the drag rate is caused 
mainly by the acoustic mode, which is lower in energy and hence easier to excite ther- 
mally, while at higher temperatures both modes contribute. At even higher temperatures 
(T > 0.67V) the separation into plasmon and non-plasmon contributions is not well-defined 
because the plasmons are heavily Landau damped and merge into the single-particle excita- 
tion continuum; cf., Fig. |2|. The Landau damping weakens the plasmon enhancement effect, 
causing the scaled r^ 1 to peak at around T « 0.57V. At extremely high temperatures, 
Tp 1 ~ T~ 3 / 2 , as shown in Appendix |C[ 

Fig. [| also shows the result of previously used expressions for the screened interaction. 
The dashed line represents the approximation that the temperature is set to zero in the 
polarization function used by Zheng and MacDonald.ll! As discussed earlier, this approxi- 
mation misses theolasmon contribution, as does the static screened approximation used by 
Jauho and Smith,H3 given by the dotted curve. Neither of these curves shows the plasmon 
enhancement demonstrated by the full numerical evaluation of Eq. Q20|). 

We have also evaluated of r^ 1 for densities which are smaller and larger than the one 
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used in the experiment of Ref. |8|. Calculations relevant for that particular density are shown 
in Ref. 14 . In Fig. [5] we show the scaled drag rate as a function of the temperature. In both 
plots we have used a well separation of 800 A. The full curves are for matched densities, 
while for the dashed (dotted) the density of layer 2 is twice (half) that of layer 1. In Fig. 
|5](a) the layer 1 density is 5. x 10 11 cm -2 , and in (b) it is 1. x 10 11 cm -2 . Since the Fermi 
temperature is thus 5 times higher in (a), the peak due to the plasmon enhancement occurs 
at a much higher temperature. Also note that the large peaks in the scaled drag rate occur 
only in the n^jnx = 1 curves. 

The insets show the drag rate as a function of the relative density ratio, which show 
clear peaks at matched densities. The peaks are signatures of the plasmon enhancement 
effect, and they occur when the Fermi velocities Vp^ of both subsystems are equal.0 This 
is because the size of the plasmon contribution to r^ 1 is determined by the smaller of the 
lm[xi{q, ^±{q))], which is roughly given by the distance between the particle-hole continuum 
and the plasmon dispersion line on the u — q plane. The min{Im[xj]} is maximized at 
matched vp^. As the vp/s start to differ from each other, the particle-hole continuum 
of the subsystem with the smaller vp, say subsystem 2, starts to move away from the 
plasmon dispersion lines (which lie above the particle-hole continuum of the subsystem 1 
with the larger Vp, 1). Consequently, the Im[x 2 (9) k>± (<?))] of the subsystem 2 is decreased 
with increasing mismatch ofvp, which accordingly decreases the total plasmon contribution 
to the drag rate. The peak at equal Fermi velocities (and hence matched densities when the 
masses of the wells are equal) should be the feature of the predicted plasmon enhancement 
which is easiest to experimentally verify. Other signatures include the aforementioned d~ 3 
behavior and the upturn in the temperature scan at approximately O.2Tp.0 



V. DRAG RATE AS A FUNCTION OF DOPANT DISTANCE FROM QUANTUM 

WELL 

Thus far, we have assumed that the r(k) is constant (which is equivalent to assuming 
that the distribution function is a shifted Fermi-Dirac function) and hence Eq. (|20|) is valid. 
This is only true if the sum of the impurity and phonon scattering rates is relatively energy 
independent and/or these rates are much smaller than the intralayer electron-electron scat- 
tering rates, as is the case when the dopants are placed approximately 700 A away from the 
quantum well. We find, however, than when the impurities are moved closer into the wells, 
Eq. ( |2ll| ) is no longer valid; in fact, when the impurities are placed 100 A from the side of 
the well, it underestimates the drag rate by approximately a factor of 2 when compared to 
Eq. flT9D because it neglects the energy dependence of the momentum relaxation time r(fc). 

To evaluate the drag rate in this case, one must be able to calculate the momentum 
relaxation rate r(fc) in the presence of impurity, phonon and electron-electron scattering. 
We adapt the formalism used to study this in 3-dimensionsIll to 2- dimensions. H We assume 
that the impurities are charged and are in an uncorrelated 5-doped layer a distance s from 
the side of the quantum well. We calculate r(k) including effects of the dynamically screened 
intralayer electron-electron, screened remote ionized impurity and screened acoustic (both 
deformation potential and piezoelectric) phonon@ scattering. Our calculations indicate that 
the r{k) is exceedingly flat when s = 700A, but as the impurities are moved in r(k) starts 
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to show a significant positive slope at around s = 300 A. This is because the charged im- 
purity potential is long-ranged, implying that it falls off rapidly with increasing momentum 
transfer,^ and hence the momentum relaxation time increases with increasing momentum. 
At large s, the electron-electron scattering, which drives the distribution to a shifted Fermi- 
Dirac function, suppresses this positive slope. However, as the impurities are moved closer to 
the well, the impurity scattering and electron-electron scattering rates become comparable, 
and hence the positive slope in r(k) develops. 

A positive slope in r(k) enhances the drag rate because it implies that there are relatively 
more carriers at higher energies (see Eq. (|12|)). These higher energy carriers are able to 
carry more current, which consequently causes increase in the drag rate. The enhancement 
is even more pronounced when we include dynamical screening because it is the high-energy 
particles which participate in the plasmon-mediated interlayer Coulomb interaction. In Fig. 
|6] we show the drag rate as a function of the charged impurities from the side of the wells 
s (which we assume is the same for both wells), for several different temperatures, all other 
parameters being fixed. The figure shows that one can obtain an enhancement in the drag 
rate of approximately a factor of 2 by moving the impurities to within 100 A of the side of 
the well, all this coming from the momentum dependence of r(k). From the figure, one can 
see the effect of the plasmons again on r^, 1 in the fact that the enhancement at small s is 
largest at T = 30K where the plasmon enhancement is the greatest, because a larger number 
of high-energy particles can now partake in the plasmon-enhanced interlayer scattering. 

Finally, since we have studied the drag rate in the case where the absolute magnitude of 
the impurity scattering rate is large (because s is small) one could ask if the diffusive form 
of the polarizability function0S0 at small q and u has any effect on our calculations at 
experimentally relevant temperatures. We conclude it does not because the crossover tem- 
perature is below which the diffusive effect should be seen goes asEl ~ 0.1Kexp[— 0.9(£/d) 2 ]. 
Since d is on the order of 400 A, and I, the mean free path, is approximately 0.1 to 1 /xm 
for s = lOOA, the crossover temperature is extremely small and experimentally irrelevant. 



VI. DISCUSSION AND SUMMARY 

Using the Boltzmann equation, we have derived a formula for the transresistivity of 
coupled two-dimensional electron gases, valid for arbitrary intralayer scattering, to lowest 
nonvanishing order in the interlayer interaction. While the Boltzmann equation does not 
include quantum effects such as weak-localization, it has been shown00 that these effects 
for systems where kpl ^> 1 are negligible. The transresistivity depends on the momentum 
relaxation times r»(fe) of both layers, which are given by the linear-response solution of the 
distribution function to an applied electric field). 

For isotropic parabolic bands, one can define a drag rate which is proportional to the 
transresistivity. In this case, in the limit where r(k) is a constant, one regains previously 
the obtained result, Eq. from our result, Eq. fll^p. A constant r(k) is equivalent to 

a shifted Fermi-Dirac distribution function under application of a small electric field. In 
experiments done previously, where the dopant layer in the modulation-doped structures 
is relatively far away from the quantum wells, r(k) is in fact relatively constant over the 
relevant energy range (i.e., within a few k^T from the chemical potential), because the 
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electron-electron scattering rate, which tends to drive the distribution to a shifted Fermi- 
Dirac function, dominates over all other scattering rates. Thus, for calculations involving 
these structures, we have used Eq. (p0|). 

The interlayer coupling is given by the screened Coulomb interaction. Due to the pres- 
ence of plasmons, both acoustic and optic, the interaction can be significantly enhanced, as 
these plasmons serve to "anti-screen" the interaction. At low temperatures, the plasmons 
are "frozen out," and hence they play no role. At higher temperatures, however, the plas- 
mons can enhance the drag rate by almost an order of magnitude over interactions which 
exclude dynamic effects of the coupled electron gas system. We have shown that the max- 
imum plasmon enhancement occurs around 30 K for GaAs when both wells are doped at 
1.5xl0 11 cm -2 . Furthermore, because the plasmons enhancement is most effective when the 
Fermi velocities of the two electron gases are equal, one should see a peak in the r^ 1 as the 
ratio of densities of the electron gases n 2 /ni is varied through 1. 

We have also investigated the effect of moving impurities closer to the quantum wells, 
which gives r(k) a positive slope. This increases the number of high-energy (and hence large 
current- carrying) particles and hence also enhances the drag rate. The effect is largest when 
the the plasmons effect is the greatest. One obtains an enhancement of approximately a 
factor of 2 when impurities are moved in from s = 400 A to s = 100 A. 
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APPENDIX A: BOLTZMAN EQUATION FOR SINGLE LAYER 



In this Appendix, we describe the formalism for the Boltzmann equation in a single 
quantum well, which we use in the main text. 



1. Formalism for intralayer scattering 



The linearized form of the Boltzmann equation for (positively charged) carriers with an 
electric field E is 

e E- (?L.) = eE-v(^f-) = -H[^{k)\. (Al) 



dp i ^ de 

From the principle of detailed balance, one can show that H is a Hermitian operator^!; i.e., 
for arbitrary functions a(k) and b(k) 

J dk a{k) H[b]{k) = Jdk b{k) H[a](k). (A2) 

Therefore, the current in the /3-direction is given by Eq. (|j) and (|A1|), 
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/dk 



-2ek B T 



(2 



dk - d Mk) (£) m 



2k B Te 2 J2E a f 

a 



dk 



v p {k) jH v a ( 



de 



(k) 



(A3) 



where the factor of 2 is for spin, and /if" is mobility tensor of the system. 



2. Impurity scattering 

We now examine a concrete example, an isotropic system dominated by impurity scat- 
tering, in which case the linearized collision operator is given by 

- H[^]{k) = k B T(^f-) (A;) r-\k) ^(fc), (A4) 



de 

where r(k) is the transport lifetime.!^ The inverse of H is trivial, 



= - iJwWY (A5) 



and hence the solution of Eq. (|A1[) , for E in the x-direction, is 



m = (A6) 



APPENDIX B: NUMERICAL EVALUATION OF THE FINITE TEMPERATURE 

POLARIZATION FUNCTION 

Here we describe an efficient technique for calculating the finite-temperature 2D-RPA-% 
for an isotropic parabolic band. The RPA x is given by 

2m r dk /°(fc + g/2)-/°(fc-q/2) 
x{q ^ ) = Wl{2,Y k.q-mu J /h- l 0+ ' (B1) 

where /° is the Fermi distribution. We define the following nondimensional quantities 

Q = q/k F , K = k/k F , Q = hu/Ep, t = k B T/E F , fr = fi/E F , x = 77^ , (B2) 

mj [Tin ) 

where k F = (27m) 1 / 2 (n is the density), and E F = h 2 k F /(2m). Note that x 1S normalized 
by the value of x(<l —> 0, uj = 0; T = 0). In these units, 
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2w J x y K x Q-Q/2-iO+ 



where 



f(K x ,K y ) = ,< K2 ^ K2 ~i /-/-\ I i ' £ = * Me^-l). (B4) 



1. Imaginary part of x 

The imaginary part of \ is given by 

Im[x(Q, = — I" dK v f°(— + 9. K v ) - f°(— --,K V ) 
1 " 2QJ-00 yJ y 2Q 2' y! J y 2Q 2' y> 



■r-i/2{—) ~ f-i/2{—. 



^-OH)'- < B5 » 



2Q 
where 

l f°° y^ 1 ^ 2 

T-x j2 {x) = —=\ dy — — , (B6) 

V^r Jo exp(y — x) + 1 

is the Fermi function of order —1/2.0 J 7 -i/2 is a well-known function, and many excellent 
approximation schemes are available to evaluate it.0 



2. Real part of x 

We use the expression by Maldagueil 

X (q, W ; T) = j Q d»> X ( Q , u; T = 0) ^—^ _ • (B7) 

Letting a± — |(g ± Q) and using the known form for u;; T = 0), we obtain 

Re[\-(Q,Si)] = -f*' 4tc03h ,| ( ;_ A , )/(2t) i [l - =^t) ^««4 - + 

" ~eM-m + 1 + Q M ' (a+) " Q M ' (< '- ) (B8) 

where 

^^r^' ^hW^MTM ' 1 "" '' 2 - (B9) 
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For any given temperature, this function M t (x) need only be evaluated once, and then a 
polynomial, Pade or a spline fit can be made to it. Then, Re[x(q, to)] can be evaluated with 
this fitted function, using Eq. ((B8|). Caution, however, should be exercised when Q <C 1 and 
Q/Q is not small. Then, the terms ±Q~ 1 M t (a±) in Eq. ( p8|) are large and approximately 
cancel, leading to large numerical uncertainly. We got around this problem by using the 
well-known asymptotic form Re[%(Q, fl)] = 2Q 2 /Q 2 , which can be derived from Eq. (|B8| ) . 
This form should be used in the regions when numerical inaccuracies plague the form given 
in Eq. flES]). 



APPENDIX C: THE LARGE TEMPERATURE LIMIT T > T F 

In the case of very large temperature, the Coulomb interactions are basically unscreened 
because the typical wavevector transfers is much larger than the debye screening wavenum- 
ber, = 27re 2 n/(e /csT), and one can assume that the interlayer interaction is the bare 
unscreened interaction. Then, under further assumptions that the carriers are nondegen- 
erate and the intralayer momentum scattering rate is constant, the we show that the drag 
rate goes as r^ 1 oc T -3 / 2 . Since r^ 1 oc T 2 for small T, the drag rate must peak at some 
temperature. 

The polarization function for a Maxwell Boltzmann distribution in the RPA is 



AW ' ; 2tq 



K 2q 2 ! K 2q 2' 



(CI) 



where Z is the plasma dispersion function.il Here as before, x = x/[ m / (F >2 ' 7T )] an d t = T/Ep, 
whereas the other variables with tildes are defined as being normalised with respect to 
temperature, 



q = hq/ \j2mk B T } 



UJ 



k B T' 



d = d^2mk B Th- 1 . (C2) 

The classical limit (q -C 1) of the polarization function Xci [which is valid here because the 
q integration is cut off at small q; see below] only depends on the ratio u/q, 

X cl {u/q) = ±-Z'(% z ). (C3) 



2t \2q 



The imaginary part is given byH 



UJ\/1X l uj 2 



Im[x d (a)/q)] = exp ( -— . (C4) 



2tq V 

From Eq. (|20|) , the drag rate for identical layers is 
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hq% ( 27re 2 m \ 2 | f°° , , 
1 ] { 1 1 dq q 6 



8n 2 nm \e irh 2 q F 

' %~ 2 sinh 2 (^/2) r^J- (C5) 

Having nonidentical layers simply leads to a change in the prefactor. Defining f D 1 as r D 1 di- 
vided by the non-temperature-dependent prefactor in the curly braces in the above equation 
yields 



Since c? oc T 1 / 2 , for large T the integrals are cut off at small q. The presence of the 
gaussian term implies in the u integral implies that this integral is also cut off at small 
values, and hence an expansion sinh 2 (w/2) ^Cj 2 /A can be made. This gives 



To' 



poo POO 

nr l / dqq~ l / dCu exp(-cu 2 /(2<f)), 
Jo Jo 

dq exp(—2qd), 

2) u (C7) 

Since t oc T and d oc T 1 / 2 , this shows that the temperature dependence of the drag rate in 
the large temperature limit, given an energy-independent intralayer momentum relaxation 
rate, is 



71 71 

7 V 2 Jo 

7TX3/2 J 



T D X OC 



In actual fact, this asymptotic behavior is reached very slowly. Our numerical evaluation 
of the drag rate indicates that for the experimental parameters of Ref. || this behavior 
of r^, 1 occurs only when T 3> lOTp, and hence is not experimentally observable for these 
parameters. 
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FIGURES 



FIG. 1. Plasmon dispersions at zero temperature for two identical wells based on Eq. fl28a|) . 
The dispersions are shown for two different values of the layer separation kpd = 4 (solid lines) and 
kpd = 12 (dashed lines), and the hatched area is the particle- hole continuum. The lower- (upper-) 
lying branch correspond to charge density oscillations in the two layers being out of (in) phase. 
The modes are well-defined only at small q. At larger q, they disappear either by merging together 
or due to Landau damping. 

FIG. 2. Scattering rates [integrand of Eq. (^C|)1 as a function of energy transfer fiuj, for 
GaAs coupled quantum wells with equal electron densities of n = 1.5 x 10 11 cm -2 , well separation 
d = Skp l = 800 A, q = 0.3 k F , and zero well widths, at two different temperatures T = 0.5 Tp and 
0.8 Tp, where Tp = 61 K. The dashed lines are the plasmon-pole approximation to the plasmon 
peaks, which is used for calculation the pole contributions in Fig. Note how the poles merge and 
are Landau damped at higher temperatures, implying that the plasmon-pole approximation is not 
valid for T > 0.6Tp and higher. 

FIG. 3. Contour of the integrand of Eq. (|20| ) for the case of two identical GaAs quantum wells 
with densities n = 1.5 x 10 n cm~ 2 , and a well separation of of d = 3.75 kp 1 = 375 A, well widths of 
200 A, as in the experiment by Gramila et all, at T = 18 K (=0.3 Tp). Adjacent contours differ 
by W~ 4 Ep/h = 8 x 10 8 s -1 , and the rightmost contour line equal to W~ 4 Ep/h. The dashed lines 
are the plasmon-pole approximation dispersion curves. 

FIG. 4. Temperature dependence of the drag rate scaled by T 2 , for the same parameters as 
in Fig. |3[ The full bold curve corresponds to calculations using the finite-T form of x(<2S w )> the 
dotted curve to using the T = form of xJ^ an d the short-dashed is based on the static screening 
approximation .111 Also shown are the plasmon-pole approximation estimates for the acoustic plas- 
mon (ap) and optic plasmon (op) contributions to the Tq , and the sum of the two (op+ap). For 
T > 0.6Tp, this approximation becomes less reliable due to large Landau damping of the modes, 
and hence we have plotted the results with dashed lines. 

FIG. 5. The drag rate as a function of temperature, for well separation d = 800 A and well 
width 200 A, for densities (a) n\ = 5 x 10 n cm~ 2 and (b) ri\ = 1 X 10 11 cm -2 . The full, dashed and 
dotted lines are for niln\ = 1, 2 and 0.5, respectively. Insets: drag rate as a function of relative 
density n^lnx, when n\ is fixed and ri2 is allowed to vary at (a) 80 K and (b) 20 K. The peaks at 
n 2/ n i = 1 are evidence of plasmon enhanced drag. 

FIG. 6. The drag rate as a function of the distance of the charged dopant layer, s, from the 
side of the well. The solid (dashed) lines were calculated using full dynamic (static) screening, and 
T = 15, 20 and 30 K. The parameters used were n = 1.5 x 10 n cm~ 2 , d = 375A, and well width 
of 100 Ain GaAs. For s ^ 400 A, the distribution functions are very close to shifted Fermi-Dirac 
functions. As s decreases, the distribution functions deviate from shifted Fermi-Dirac functions, 
and the drag rate increases markedly. 
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